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The binary Bi-chalchogenidcs, Bi2Cli3, are widely re- 
garded as model examples of a recently discovered new 
form of quantum matter, the three-dimensional topolog- 
ical insulator (TI) [1-4]. These compounds host a sin- 
gle spin-helical surface state which is guaranteed to be 
metallic due to time reversal symmetry, and should be 
ideal materials with which to realize spintronic and quan- 
tum computing applications of TIs [5] . However, the vast 
majority of such compounds synthesized to date are not 
insulators at all, but rather have detrimental metallic 
bulk conductivity [2, 3]. This is generally accepted to re- 
sult from unintentional doping by defects, although the 
nature of the defects responsible across different com- 
pounds, as well as strategies to minimize their detri- 
mental role, are surprisingly poorly understood. Here, 
we present a comprehensive survey of the defect land- 
scape of Bi-chalchogenide TIs from first-principles calcu- 
lations. We find that fundamental differences in the ener- 
getics of native defect formation in Te- and Se-containing 
TIs enables precise control of the conductivity across the 
ternary Bi-Te-Se alloy system. From a systematic angle- 
resolved photoemission (ARPES) investigation of such 
ternary alloys, combined with bulk transport measure- 
ments, we demonstrate that this method can be utilized 
to achieve true topological insulators, with only a single 



Dirac cone surface state intersecting the chemical poten- 
tial. Our microscopic calculations reveal the key role of 
anti-site defects for achieving this, and predict optimal 
growth conditions to realize maximally-resistive ternary 
TIs. 

Figure 1 shows the calculated Fermi-level-dependent 
formation energies of native defects in the binary TIs 
Bi 2 Se3 and Bi 2 Te3. For Bi 2 Se3 under Bi-rich/Se-poor 
conditions (Fig. 1(a)), Se vacancies in the two inequiv- 
alent chalcogen layers, Vg c and Vg e , are the dominant 
donor defects. In particular, Vg e has the lowest forma- 
tion energy of all of the native defects for Fermi levels 
across the entire bulk band gap, and even up in to the 
conduction band. This will result in a strong propensity 
for the formation of n-type defects, without significant 
compensation from p-type defects for bulk Fermi levels 
up to at least 0.1 eV above the conduction band minimum 
(CBM). This is entirely consistent with our experimental 
measurements shown in Fig. 2. Our ARPES measure- 
ments (Fig. 2(a)), which probe the occupied electronic 
structure, show not only the bulk valence bands and the 
topological surface state, but also occupied bulk conduc- 
tion band states. The Fermi level is located a little over 
0.1 eV above the conduction band minimum (CBM), as 
expected from our calculations, resulting in a large n- 
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FIG. 1. Formation energies, as a function of Fermi level rela- 
tive to the VBM, of donor (solid lines), acceptor (dashed lines) 
and electrically-inactive (dot-dashed lines) defects in Bi2Se3 
under (a) Bi-rich and (b) Bi-poor conditions. (c,d) Equiv- 
alent calculations for Bi2Te3. (e) Calculated valence band 
offset and resulting band alignment of Bi2Se3 and Bi2Te3. 



type conductivity, and a temperature-dependent resis- 
tivity (Fig. 2(b)) characteristic of a metal rather than an 
insulator. 

For comparison, we also calculate the formation ener- 
gies for native defects in Bi2Se3 under Bi-poor/Se-rich 
conditions (Fig. 1(b)). Compared to the defect energet- 
ics under Bi-rich conditions, one would naively expect 
the formation energy of Vs e to increase, while that of the 
acceptor- type Bi- vacancy to decrease, and this is indeed 
seen in our calculations. On this basis alone, much lower 
residual rt-type conductivities could be expected under 
Bi-poor growth conditions, with the Fermi level moving 
into the bulk band gap (tending towards the intersec- 
tions of the formation energies of Vg e and Vsi)- We find, 
however, that the donor-type Se anti-site defect, Seei, 
becomes the lowest energy defect throughout the band 
gap. Therefore, as well as the commonly assumed Vs c , 



TABLE I. Anti-site disorder energies, A_E a( j = 
| [£a b + E Ba - 2E pulc ], where E Ab is the total energy 
of a supercell containing an Ab defect, and E puTC is the 
energy of the stoichiometric Bi2Cli2Ch 2 supercell. 
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our calculations indicate that SeBi can play a significant 
role in driving the unintentional conductivity of Bi2Se3. 
For all possible growth conditions, the lowest energy de- 
fect is a donor (either Vg e or Seei), even when the Fermi 
level lies at, or slightly above, the CBM. This explains 
why crystals of Bi 2 Se3 always display unintentional n- 
type conductivity, which can only be compensated by 
suitable extrinsic p-type doping [6]. 

The defect physics of Bi2Te3 is rather different, 
with anti-sites being the dominant defects under 
both Bi-rich/Te-poor and Bi-poor/Te-rich conditions 
(Fig. l(c,d)). Under Bi-poor conditions (Fig. 1(d)), the 
Teei donor defect has the lowest formation energy of all 
of the native defects. As for Bi2Se3 under Bi-poor con- 
ditions, this will yield unintentional n-typc conductivity 
as often observed in experiment [3] , with significant com- 
pensation by acceptor Vbi centres not expected for Fermi 
levels within the bulk band gap. However, in contrast 
to Bi2Se3, under Bi-rich conditions, the formation en- 
ergy of the acceptor Bi anti-site defect, Bi^ e , becomes 
smaller than that of the chalcogen vacancies (Fig. 1(c)). 
This promotes a natural tendency for unintentional p- 
type conduction when Bi2Te3 is grown under Bi-rich con- 
ditions, consistent with both single- crystal growth exper- 
iments [7] as well as recent studies on MBE-grown thin 
films [8]. Thus, our calculations reveal that the defect 
landscape of Bi2Te3 is dominated by anti-site defects for 
all growth conditions, and anion vacancies play a much 
less significant role than in Bi2Se3. 

In fact, as shown in Table I, the energy cost of anti-site 
disorder in Bi2Te3 is approximately half that of Bi2Se3. 
This is due to the more similar ionic radii of Bi and Te [9] , 
and the relatively small differences in electronegativity 
between the two species: the resulting anti-site defects 
represent a relatively low-energy configuration. We note 
that these antsites serve to oppose the "expected" po- 
larity of the materials, yielding p-type conduction under 
typical n-type growth conditions, and n-type samples un- 
der typical p-type conditions. While growth conditions 
between these two extremes can yield defect energetics 
which drive the Fermi level into the bulk band gap, we 
note that the small size of this energy gap will make 
achieving robust insulating behaviour very difficult for 
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FIG. 2. (a) ARPES measurements of 
layered Bi-chalchogonides. The bulk 
valence bands (BVBs) and topologi- 
cal surface state (TSS) are clearly ob- 
served for all compounds. Occupied 
bulk conduction band (BCB) states are 
only observed for Bi2Se3 and Bi2Se2Te, 
indicating a transition from degener- 
ately doped semiconductors to true 
bulk TIs as the Te content is in- 
creased. Temperature-dependent resis- 
tivity measurements (b) confirm this 
trend, showing metallic and insulating 
behaviour for Bi2Se3 and Bi2Te2Se, re- 
spectively. The inset shows a schematic 
representation of their electronic struc- 
ture, with the Fermi level located in the 
conduction band for Bi2Se3 and within 
the bulk band gap for Bi2Te2Se. 



this compound. 

The qualitative differences in the native defect be- 
haviour of Bi2Se3 and Bi2Te3 can be understood from 
their band alignment. We have computed the natural 
valence band offsets of these materials using the method- 
ology of Zunger and co workers [10, 11]. We find a stag- 
gered "type II" [12] offset (Fig. 1(e)), with the valence 
band maximum (VBM) of Bi2Te3 0.26 eV higher in en- 
ergy than that of Bi2Se3. The smaller ionization poten- 
tial of Bi2Te3 suggests an increased preference for hole 
formation [13], which fully supports our microscopic cal- 
culations. The CBM of Bi2Se3, on the other hand, is 
only O.lOeV below that of Bi2Te3, explaining why both 
materials display similar n-type behaviour under Bi-poor 
growth conditions. 

Within a conventional semiconductor band engineer- 
ing methodology, this suggests that alloying Bi 2 Se 3 and 
Bi2Te3 could be a suitable way to realize bulk insula- 
tors, where the topological surface state conduction is no 
longer shunted by a large residual bulk conductivity. In- 
deed, transport measurements of the ternary compound 
Bi2Te2Se have already found a much more insulating bulk 
resistivity than for the binary compounds [14], although 
previous ARPES measurements still showed the occupa- 
tion of a small number of states at the bottom of the 
conduction band [15]. In contrast, our APRES measure- 



ments (Fig. 2(a)) show that, upon moving towards Te- 
rich Bi-Te-Se alloys, the conduction band is readily de- 
pleted of carriers. In particular, for both Bi 2 Te2Se and 
Bi2Te2.5Seo.5, only the topological surface state inter- 
sects the chemical potential, as desired for a true TI [16]. 
Both of these compounds exhibit a temperature depen- 
dence of their resistivity indicative of bulk insulators. For 
Bi 2 Te 2 Se (Fig. 2(b)), the low-temperature resistivity is as 
much as two-to-three orders of magnitude higher than in 
Bi2Se3. However, Bi2Te2.5Seo.5 is approximately a factor 
of 6 less resistive than Bi 2 Te2Se, due to its smaller band 
gap. 

Thus Bi2Tc2Se can be seen as a more ideal TI, and we 
perform explicit calculations for this compound in order 
to elucidate the microscopic origin of its enhanced resis- 
tivity as compared to the binary compounds. Both the Te 
and Se chemical potentials can be simultaneously varied, 
subject to the constraints that Bi2Te2Se has lower forma- 
tion enthalpy than binary compounds of Bi and Te/Se, 
or than elemental Bi, Te, or Se. Considering these limits, 
we follow the approach of Walsh et al. [18] and Persson 
et al. [19] to calculate a phase diagram for growth of 
Bi 2 Te2Se, shown in Fig. 3(a). We consider five repre- 
sentative environments in which to calculate formation 
energies of native defects, shown in Fig. 3(b). 

In addition to the defects considered for the binary 
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FIG. 3. (a) Calculated phase diagram of Bi2Te2Se as a function of Te and Se chemical potential, (b) Formation energies for 
defects in Bi2Te2Se, computed for the five sets of chemical potentials shown in (a). 



compounds, it is now possible to have anion-on-anion 
anti-site defects (Sex e and Tes e )- In fact, we find that 
these centres have the lowest formation energy of all na- 
tive defects across the entire phase diagram. While they 
have no transition levels within the band gap, and so 
are electrically inactive, the energy barrier for anion-on- 
anion anti-site disorder is as low as only 0.13 eV (Table I). 
It therefore seems inevitable that there will be a signif- 
icant level of anti-site disorder in these compounds, and 
Bi2Te2Se alloys will not form as ordered structures as is 
commonly assumed. 

Depending on the growth environment, the lowest en- 
ergy p-type defects are the Bi vacancy or the Bix e anti- 
site, with the dominant n-type defects being the Te va- 
cancy or the Teei anti-site. Se vacancies and anion-on- 
cation anti-sites invariably have higher formation energy 
than their Te counterparts. In all cases, the formation 
energy of the dominant donor and acceptor defects cross 
over within the band gap. Under equilibrium conditions, 
the Fermi level will tend to be pinned close to this cross- 
ing point, represented by the vertical red dotted lines 
in Fig. 3(b). For Te/Se-rich conditions (C and D in 
Fig. 3(b)), this is very close to the CBM, and so the 
bulk conductivity will likely still be rather high. How- 
ever, for Te- and Se-poor growth conditions (A, B, and E 
in Fig. 3(b)), this level lies close to the middle of the band 
gap, where the formation energy for the doubly-charged 
donor Vr c crosses that of the singly-charged acceptor 
Bixo- Our calculations therefore indicate that these (par- 
ticularly point A) represent the ideal set of growth con- 
ditions in which to realize maximally-resistive Bi2Te2Se. 
Only considering vacancies, the effective Fermi level pin- 
ning would be shifted towards much more n-type condi- 
tions, and so it is again clear that anti-site defects play a 
key role controlling the unintentional bulk conductivity 



of ternary, as well as binary, topological insulators. 



EXPERIMENTAL SECTION 

First principles calculations: Density Functional The- 
ory (DFT) calculations were performed using the projec- 
tor augmented wave method [20] implemented within the 
VASP code [21]. Exchange and correlation were treated 
within the PBE functional [22], using a planewave cut- 
off of 300 eV and a fc-point sampling of 10 x 10 x 10 for 
the 5 atom tetradymite unit cell. The structure was 
deemed to be converged when the forces on all of the 
atoms were less than 0.01 eV A -1 . Defect calculations 
were performed using the method described in Ref. [4], 
using 4x4x1 expansions of the hexagonal representation 
of the unit cell (i.e., 240 atom supercells) with a 2x2x1 
Monkhorst-Pack special fc-point grid. All calculations in- 
clude spin-orbit coupling. 

Experimental details: Bi2Se3, Bi2Se2Te, Bi2Te2Se, and 
Bi2Te2.5Seo.5 crystals were prepared by melting high pu- 
rity elements (5N) of Bi, Se, and Te in the ratios 2:3:0, 
2:2:1, 2:1.05:1.95, and 2:0.5:2.5, respectively. ARPES 
measurements were performed using Scienta R4000 hemi- 
spherical analysers at beamline 5-4 of the Stanford Syn- 
chrotron Radiation Lightsource (SSRL) and the CAS- 
SIOPEE beamline of synchrotron SOLEIL. The pho- 
ton energies were between 14 and 18 eV, and the sam- 
ple temperature was ~ 10 K. Samples were cleaved at 
the measurement temperature in a pressure better than 
3 x 10~ n mbar. Further theoretical and experimental 
details are given in supplementary information. 
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SUPPLEMENTARY INFORMATION 
FIRST PRINCIPLES CALCULATIONS 

Prior to performing supcrccll calculations of the de- 
fect formation energies, we calculated the bulk electronic 
structure of stoichiometric Bi 2 Se3, Bi 2 Te3, and Bi 2 Te2Se 
using the paramaters as indicated in the manuscript. 
These calculations yielded indirect band gaps of 0.29 eV 
and 0.27 eV for Bi 2 Se3 and Bi 2 Te 2 Se, respectively, consis- 
tent with experiment. For Bi 2 Tc3, our calculated indirect 
band gap of 0.08 eV is approximately half of the exper- 
imental band gap [1-3]. We corrected this error for the 
formation energy calculations using a scissors operator. 
For Bi 2 Te 2 Se, we assume an ordered alloy in the calcu- 
lations. Namely, starting from Bi 2 Te3, this compound is 
realized by replacing the middle chalcogen layer of the 
quintuple- layer structure (Ch 2 in Fig. 4(a)) with Se. 

We calculate defect formation energies using the 
method described elsewhere [4] , incorporating corrections 
for the finite size of the supercell [5] as well as for band- 
filling resulting from shallow donors [6]. We have con- 
sidered the formation of all isolated native defects (see 
Fig. 4(b) for site definitions): bismuth vacancies (Vei), 
both possible anion vacancies (Up h and VqjJ, cation and 
anion interstitials in the tetrahedral and octahedral con- 
figurations (Bi/Chf* and Bi/Ch° ct ), both Bi on chalco- 
gen antisites (Big h and Bi^ h ), and the Ch on Bi antisite 
(Chei). For each compound, we perform our calculations 
under the limit of Bi-poor and Bi-rich conditions, respec- 
tively. Physically, this corresponds to varying the partial 
pressures during growth, which can be achieved by ad- 
justing the relative Bi and Ch flux as well as the substrate 
temperature during molecular-beam epitaxy (MBE), or 
by adjusting the Bi/Ch ratio in the starting mixture for 
bulk crystal growth. For Bi-rich/Ch-poor conditions, we 
found the bounding chemical potential values to be set 
by BiCh (rather than Bi 2 Ch 3 ) formation. 

SAMPLE GROWTH DETAILS 

The samples were synthesized starting from powder 
mixtures in the ratios given in the main manuscript. The 
mixtures were reacted in sealed, evacuated quartz tubes 
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FIG. 4. (a) Quintuple-layer (QL) building block of the I^Ch^Ch 2 (Ch = Se,Te) structure, (b) The tetrahedral and octahedral 
interstitial sites between the quintuple layers. 



at 850°C for 2 days. For Bi2Se3, this was followed by 
cooling at 2-3° C/h to 650° C and annealing at this tem- 
perature for 7 days before quenching to room tempera- 
ture. The Bi2Te2Se, Bi2Te2.sSeo.5, and Bi2Se2Te were 
cooled at 2° C/h to 450° C, followed by cooling at 50° C/h 
to room temperature. The Bi 2 Te 2 Se sample was then 
resealed in a quartz tube and annealed at 600° C for 2 
weeks to obtain the insulating behaviour shown in Fig 2 
of the main manuscript. 
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